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Abstract. 

We present optical and near-infrared observations of the afterglow of the gamma-ray burst GRB 050904. We derive a photomet- 
ric redshift z = 6.3, estimated from the presence of the Lyman break falling between the / and J filters. This is by far the most 
distant GRB known to date. Its isotropic-equivalent energy is 3.4x 10 53 erg in the rest-frame 1 10-1 100 keV energy band. Despite 
the high redshift, both the prompt and the afterglow emission are not peculiar with respect to other GRBs. We find a break in 
the ./-band light curve at f b = 2.6 ± 1.0 d (observer frame). If we assume this is the jet break, we derive a beaming-corrected 
energy E y ~ (4 -f 12) x 10 sl erg. This limit shows that GRB 050904 is consistent with the Amati and Ghirlanda relations. This 
detection is consistent with the expected number of GRBs at z > 6 and shows that GRBs are a powerful tool to study the star 
formation history up to very high redshift. 
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* Based on observations carried out at ESO telescopes under pro- 
gram Id 075.D-0787. We also used data from Telescopio Nazionale 
Galileo and Centra Astronomico Hispano Aleman. 



1. Introduction 

Gamma-ray bursts (GRBs) are intense, short pulses of gamma 
rays, occurring at random positions in the sky. They emit large 
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amounts of energy (up to ~ 10 53 erg) and thus are detectable up 
to cosmological distances, possibly at z = 10 and beyond. The 
discovery of high-redshift GRBs is one of the main goals of the 
Swift mission (Gehrels et al. 2004) and several ground-based 
facilities were prepared to be able to detect their counterparts 
(e.g. Chincarini et al. 2003). 

GRBs are accompanied by long-lasting X-ray, optical, and 
radio counterparts (afterglows). In the first few hours after the 
explosion, these are much brighter than any known quasar. 
In particular, optical and near-infrared (NIR) afterglow spec- 
troscopy can easily provide their redshifts and probe the inter- 
vening gas along the line of sight. The possible biases affecting 
the GRB distribution are likely very different from those af- 
fecting other classes of objects, such as AGNs, damped Lya 
absorbers, and galaxies. Therefore GRBs and their afterglows 
provide an effective and independent way to probe the high- 
redshift Universe. As a matter of fact, the observed redshift 
distribution of long-duration GRBs, before the Swift era, p eaks 
at z ~ 1.6 and extends up to z — 4.5 JAndersen et allEoOOh . A 
significant fraction of GRBs has been p redicted t o occur even at 
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vides a direct measurement of the star formation rate in the 
early Universe. Finally, high-redshift GRBs may also serve as 
useful tools to study the geom etry of the Universe JAmati et all 
l2002HGhirlanda et alll20ol . 

Our group, the MISTICI 1 collaboration, has been pursu- 
ing follow-up studies of GRB afterglows for several years, 
concentrating on various aspects of the GRB phenomenon, 
among them the study of their host galaxies and of foregr ound 
DLA systems at high redshift jFiore et all 120051: iD'Eliaet all 
2005]). Here we present optical and NIR photometric observa- 
tions of GRB 050904. This burst was detecte d by the Burst 
Alert Telescope (BAT; iBarthelmv etail l2005h onboard Swi ft 
on 2005 September 4 at 01:51 :44 UT JCummings et all 120051) . 
It was a lo ng, bright burst with an observed duration 7go = 
225 ± 10 s isakamoto et all 120051) . The 15-150 keV fluence 
was (5.4 + 0.2) x 10~ 6 erg cm" 2 . The BAT spectrum can be 
described by a hard power law with photon index « 1.2. Long- 
lastin g flaring activity was detected in the X-ray band up to 
50 ks JCusumano et allEoollWatson et alll2005l). A high red- 
shift va lue for this GRB was first suggested by lHaislip et all 
d2005ah. determined photometrically by us (z = 6.1 + HZ; 



Antonelli et al., 2005) and then confirmed spectroscopically by 
Kawai e t al. (2005), who measured z = 6.29 + 0.01. A photo- 



met ric redshift has bee n derived also bv lHaislip et all (|2£)05b) 
and lPrice et al. 1 ll2005t) . To date, this is by far the most distant 
GRB discovered and one of the most distant objects known in 
the Universe. 



1 Multiwavelength Italian Swift Team with International Co- 
Investigators. 



Fig. 1. Three colour-composite image based on our NIR bands 
(J, H, and K) obtained with VLT+ISAAC on 5 September. The 
circle represents the 6" XRT afterglow error box. The lines 
show the position of the optical afterglow. The field is 65"x65" 
wide. North is up, East is left. 



2. Observations and data analysis 

Our observations of the field of GRB 050904 were mostly 
performed with the 8.2m ESO Very Large Telescope (VLT), 
equipped with the ISAAC and FORS2 instruments (Table Q. 
Further data were collected with the 2.2m and 3.5m Calar Alto 
telescopes (equipped with CAFOS and LAICA, respectively), 
the 2.5m Isaac Newton Telescope (equipped with WFC), and 
the 3.6m Telescopio Nazionale Galileo (equipped with NICS). 
A NIR counterpart to GRB 050904 was first reported by Haislip 
et al. (2005a). In our first images, taken on 2005 Sep 4/5, 
this object was not present. A pointlike source was however 
clearly vis ible in the IzJHKg bands, at the edge of the XRT 
error box dMineo et all 120051) . Its coordinates were anooo = 
00 h 54 m 50 s .83, (5j2ooo = +14°05'1070. We suggested this source 
to be the aft erglow of GRB 050904 (D' Avanzo et al., 2005; see 
also Fig. 0. Nvsewander et al. I J2005h subsequently provide d 
revised coordinates for the object of lHaislip et all (12005a). 
showing it was coincident with our candidate. The afterglow 
was also detected by the TAROT robotic telescope up to eight 
minutes after the GRB (Klotz et al. 2005). 

Data reduction was carried out following standard proce- 
dures. Aperture photometry was performed using SExtractor 
(Bertin & Arnouts 1996). Flux calibration was achieved by ob- 
serving Landolt standard fields (/ band) and against the SDSS 
(z band) and 2MASS (NIR bands) surveys. 

In Fig.|2]we show our data, together with those by Haislip 
et al. (2005b). The light curve in the J band, for which we have 
the best coverage, shows a complex behaviour. Haislip et al. 
(2005b) have reported a flattening at early times, suggesting 
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Table 1. Observation log and photometry of the transient 
source. Errors are at the Icr confidence level, while upper lim- 
its are given at 3<x. Note that the CAFOS and FORS2 /-band 
filters (which we called I t and I 2 ) are different (their central 
wavelengths being 8500 A and 7680 A, respectively). Since the 
Lyman a break falls within this band dropout is occuring inside 
this band, this leads to a significant magnitude difference. 



Date 


Instr. 


t-t 


Exp. 


Filter 


Magnitude 


(UT) 




(d) 


(min) 






05/09/05 04:47 


ISAAC 


1.125 


15 


K s 


18.24+0.07 


05/09/06 05:23 


ISAAC 


2.146 


30 


K s 


19.12+0.07 


05/09/05 04:22 


ISAAC 


1.104 


20 


H 


19.09+0.07 


05/09/06 04:42 


ISAAC 


2.125 


30 


H 


19.96+0.07 


05/09/05 02:31 


NICS 


1.029 


50 


J 


19.58+0.14 


05/09/05 03:55 


ISAAC 


1.087 


20 


J 


19.92±0.04 


05/09/06 04-01 


ISAAC 


2.092 


35 


j 


20 75+0 07 


05/09/07 04:23 


ISAAC 


3.104 


60 


J 


21.66±0.08 


05/09/08 05:44 


ISAAC 


4.162 


60 


J 


21.91±0.08 


05/09/09 09:26 


ISAAC 


5.317 


60 


J 


22.45±0.19 


05/09/11 06:18 


ISAAC 


7.183 


102 


J 


> 23.2 


05/09/05 02:25 


WFC 


1.021 


28 


z 


21.50±0.25 


05/09/05 07:15 


FORS2 


1.250 


90 


z 


21.80±0.17 


05/09/04 23:17 


CAFOS 


0.892 


50 


h 


21.89+0.20 


05/09/05 01:35 


CAFOS 


0.987 


50 


h 


22.03+0.20 


05/09/05 03:38 


CAFOS 


1.075 


50 


h 


22.40±0.20 


05/09/05 06:54 


FORS2 


1.208 


60 


h 


24.10+0.18 


05/09/05 01:30 


CAFOS 


0.983 


90 


R 


> 24.1 


05/09/05 02:48 


LAICA 


1.038 


30 


V 


> 24.1 
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Fig. 2. Light curve of the afterglow of GRB 050904, together 
with the best fit (dashed lines) computed in the J band. The 
same fit is reported for the other bands. Open symbols repre- 
sent data from Haislip et al. (2005b). The vertical dotted line 
marks the epoch at which we computed the spectral energy dis- 
tribution (Fig. 3). 



that the early emission could be related to a different compo- 
nent, possibly the reverse shock. Our data extend for signifi- 
cantly longer and show a further steepening. Excluding the data 
at t < 0.4 d, we can model the whole dataset with a smoothly 
broken powerlaw: F(t) = 2F b/[(t/tb) a ' +(t/tb) a2 ], where a u a 2 
are the early- and late-time slopes, % is the break time and F b 
is the flux at t b . The fit provides a\ = 0.72+° ^, a 2 = 2.4 + 0.4, 
fb = 2.6+ 1 .0 d. Errors are at the Icr confidence level throughout 
the paper. 

Figure[3]shows the spectral energy distribution of the after- 
glow 1.155 d after the GRB (the epoch around which our mea- 
surements cluster). A marked dearth of flux is observed blue- 
ward of the J passband, while the J—H and H—K colors are not 
particularly red. Dust extinction cannot produce such a sharp 
cutoff, so that a redshifted Lyman dropout must be invoked (due 
to the strong absorption by neutral Hydrogen at wavelengths 
shorter than Lyff). The redshift of GRB 050904 was computed 
using two different programs. The first is a code based on a 
widely used technique to derive photometric redshifts of galax- 
ies and quasars (e.g. Fontana et al. 2000), modified to take into 
account the afterglow properties. In this case we computed a li- 
brary of synthetic GRB afterglow spectra at arbitrary redshifts, 
modeled with power laws (F v oc v ^). Reporting all our pho- 
tometric measurements at a common epoch (adopting the ob- 
served decay law), the best fit gives z = 6.3 + 0.07 and - 
1.25+ 0.25. The second is the z-ph-REM code ( F ernandez-Sotol 
120041) . specifically designed for GRB afterglows and based on 
the work by Fernandez-Soto et al. (1999). This code assumes 
that the temporal and spectral dependence of the afterglow flux 
follows a power law (F v oc r a v~P~). In this case, als o includ- 
ing the UVOT optical limits JCucchiara et all I2005I) . the best 
fit provided z = 6.3 + 0.1 and /3 — 1.2 + 0.3. Both codes use a 
^-minimisation technique to find the best-fitting spectral tem- 
plate to the observed colors. We note that t he spectral i ndex i s 
significantly redder than that reported by IPrice et al~l lEc 
which was however measured at an earlier epoch. 



3. Discussion 

Despite its high redshift, the optical afterglow of GRB 050904 
is not peculiar with respect to other GRBs. For example, its 
/?-band spectral luminosity extrapolated at t = 12 hr (rest- 
frame time) is 2 L v ~ 10 31 erg s~' Hz -1 , consisten t with the 
values typically observed fo r previous GRBs JNardini et all 
2005; Liang & Zhang! Eooll) . Its spectral and temporal indices 
are also typical among GRB afterglows. In particular, the after- 
glow colours leave small room for dust extinction, if any (note 
that we are observing the rest-frame ultraviolet band). The pre- 
cise gamma-ray energy budget of this burst is not known, since 
the observed peak spectral energy lies outside the sensitivity 
BAT range. Cusumano et al. (2005) infer that the isotropic- 
equivalent gamma-ray budget E y ^ was between 6.6 x 10 53 
and 3.2 x 10 54 erg, depending on the bolometric correction 
(E rM - 6 x 10 52 erg in the 20-100 keV rest-frame band). The 
steepening in the light curve may be due to a jet break. In this 



Assuming a cosmology with Q. m = 0.27, Ci A = 0.73, and ho = 
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Fig. 3. Photometric spectral energy distribution of the after- 
glow. Magnitudes were reported at a common epoch (1.155 d 
after the GRB) using the measured decay law. The solid line 
shows the best fit model for the afterglow. Data are corrected 
for Galactic extinction (Ay = 0.21 mag). 



case, adopting the standard formalism JSari. Piran & Halpernl 
Il999l) . we can infer a jet half-opening angle i?j et ~ 3°, assuming 
a radiative efficiency 77 = 20% and a circumburst medium den- 
sity n — 3 cm 4 (as in Ghirlanda et al. 2004). In this case the 
beaming-corrected energy would be E y ~ (4 4- 12) x 10 51 erg. 
For this burst the rest-frame peak energy is constrained to be 
E p k, 1100 keV. In the planes E p vs Ey; lso and E p vs E y , 
the derived limits show that GRB 05090 4 is consistent with 
both the Amati and G hirlanda relations jAmati et all [2002; 
Idhirlanda et ail 12004 . 

The optical/NIR afterglows of GRBs detected by Swift are 
dimmer than those of other missions, on the average by 2.5 mag 
at t — 12 hr (Berger et al. 2005). More than one third have no 
optical counterpart, despite having a well-localized X-ray af- 
terglow. For a good fraction of these bursts, this could be due 
either to a high redshift and/or to dust extinction (e.g. Zhang & 
Meszaros 2004 and references therein; see also Filliatre et al. 
2005). In about 10 months of operations, the number of Swift 
long-duration GRBs with measured redshift has been measured 
is 16, while the number of non-Swift GRBs with redshift is 43 
(obtained over a time period of 7 years). The redshift distribu- 
tion of the pie-Swift bursts has a peak at z ~ 1 .6, extending to 
higher values (the record was z = 4.50; Andersen et al. 2000). 
Although the number is still small, the 16 Swift long-duration 
GRBs with measured redshift have a statistically higher red- 
shift, with an average value of 2.8 (with one object at z > 6, 
two with 4 < z < 6, three with 3 < z < 4 and four with 
2 < z < 3, e.g. Jakobsson et al. 2005). It is very likely that 
some of the Swift bursts without optical counterpart are also 
at high redshift. Moreover, the number of bursts for which it 



is possible to perform low- and high-resolution spectroscopic 
studies at moderate redshift (z > 2) is rapidly increasing (e.g. 
D'Elia et al. 2005; Chen et al. 2005). 

Finally we note that a GRB survey is much more effi- 
cent in finding high-redshift objects than those targetting other 
sources. For instance only seven quasars out of ~ 76 000 (se- 
lected from a sample of ~ 850 000 sources for which spectra 
were acquired, over a 4783 deg 2 region of the sky) were de- 
tected in the Sloan Digital Sky Survey 3 at z > 6. For compari- 
son, Swift is monitoring about 1 .4 sr of the sky at any time, and 
it has been able to catch one event at z > 6, out of a sample 
of 16 with measured redshift, selected from a sample of ~ 60 
GRBs. 

4. Conclusions 

We observed the afterglow of GRB 050904 in the VRIzJHK, 
filters. We derived a photometric redshift z = 6.3 + 0.1 using 
two independent techniques. This redshift has been estimated 
by the Lyman break falling between the / and J filters, therefore 
it is very robust an d has been confirmed spectroscopically by 
iKawai et al] J2005). who measured z = 6.29 + 0.01. 

The high redshift of GRB 050904 opens the question of 
whether its progenitor was a Pop III star (see, for example, 
Woosley & Heger 2004; Scannapieco et al. 2005). In this case, 
however, it would be difficult for the star to lose its envelope. 
Wind-driven mass loss depends upon metallicity, and it scales 
asMcc z - 5 * - 9 (e.g. Nugis & Lamers 2000; Kudritzki 2002; 
Vink & de Koter 2005). The discovery of GRB 050904, possi- 
bly in a low metallicity environment like the one characterizing 
the early Universe, hints at the need for a close companion that 
can efficiently strip off the envelope of the progenitor. However, 
the efficiency of producing GRBs by Pop III stars is currently 
unknown and may be very low (e.g. Bromm & Loeb 2005). 

A more convincing alternative is the possibility that the pro- 
genitor of GRB 050904 was a massive, non-pristine star. This 
implies that star formation was already active at z ^ 6, and that 
metal enrichment had already started. This finding is consistent 
with recent measurement s of the comoving lumi nosity density 
of star-forming galaxies ( Giavali sco et alll2004 . which show 
that the star formation rate var ies only slowly with re dshift over 
the range 2 < z < 6 (see also lLe Fevre et alll2005l) . These re- 
sults are nevertheless still subject to important - and not well 
known - corrections due to incompleteness, dust absorption, 
and systematic effects. GRBs are not affected by dust absorp- 
tion, and are thus in principle very effective in tracing star for- 
mation even in this high redshift range. 

Swift has been able to catch one GRB at z ~ 6, out of a 
sample of 16 events for which the redshifts have been mea- 
sured. By using the whole sample of about 60 GRBs discovered 
by Swift, we find, as a simple application of Poisson statistics 
in the small-number regime (Gehrels 1986), that 4*? events 
are expected to be discovered at z > 6. A fraction of about 
4/60 ~ 7%, as we have derived on empirical grounds, is consis- 
tent with the theoretical estimate recently provided by Bromm 
& Loeb (2005), who predicted that about 10% of the events dis- 
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covered by Swift are at z > 6. Once we will be able to measure 
a significant number of high GRB redshifts, we will be in a po- 
sition to sample the star formation rate of the Universe within 
the first Gyr from the Big Bang. 
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